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Abstract 

Electron transport through parallel double quantum dot system with interdot tunneling and 
strong on-site Coulomb interaction is studied in the Kondo regime by using the fmite-C/ slave boson 
technique. For a system of quantum dots with degenerate energy levels, the linear conductance 
reaches the unitary limit (2e 2 /h) due to the Kondo effect at low temperature when the interdot 
tunneling is absent. As the interdot tunneling amplitude increases, the conductance decreases in 
the singly occupied regime and a conductance plateau structure appears. In the crossover to the 
doubly occupied regime, the conductance increases to reach the maximum value of G = 2e 2 /h. 
For parallel double dots with different energy levels, we show that the interference effect plays an 
important role in the electron transport. The linear conductance is shown to have an asymmetric 
line shape of the Fano resonance as a function of gate voltage. 

PACS numbers: 72.15.Qm, 73.21-La, 73.40. Gk 
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I. INTRODUCTION 



matter systems, including an impurity atom on metal surface jj], single-electron 
5, 6], quantum dot in AB interferometer 0, 0]. 



Due to the wave nature of electrons and the confined geometries in mesoscopic sys- 
tems, the interplay between interference and interaction becomes one of the central issues 
in mesoscopic physics. Preservation of quantum coherence in electron transport through 
an interacting regime has been manifested in the observed Kondo effect in semiconductor 
quantum dot systems|l|, and more explicitly in the conductance AB oscillation in the inter- 
ference experiment with a quantum dot embedded in one arm of an Aharonov-Bohm(AB) 
interferometer P). Recently, the Fano resonance has attracted much research interest as 
another important interference effect in mesoscopic systems . The Fano effect was first pro- 
posed as a result of the interference between resonant and non-resonant processes in the field 
r"i 

of atomic physics |3|]. It is found to be a ubiquitous phenomenon observed in a large vari- 
ety of experiments including neutron scattering, atomic photoionization, Raman scattering, 
and optical absorption. One recent progress is the observation of the Fano resonances in 
condensec 

In this paper we show the Fano effect, which can be manifested by gate voltage depen- 
dence of the linear conductance, is also important for the electron transport through double 
quantum dots(DQDs) in parallel configuration. For electron tunneling through quantum 
dots, it is well known that the strong on-site coulomb interaction leads to the Kondo effect 
at low temperatures, so that the coexistence of the Fano resonance with the Kondo effect is 
expected to yield interesting transport phenomena. Electron transport through DQDs with 
series q| and parallel^] configurations have been realized in experiments, through which 
studies on the molecular states of the double dots and also the interference effect are carried 
out. Most of theoretical studies ll[ [3, [kJ [lj] are devoted to electron transport through 
DQDs connected in series, while relatively little attention is paid for the parallel configu- 
ration case especially for the system in the Kondo regime Q]. For the DQD system with 
a parallel coupling, interference effect should play an important role. Thus, in order to 
understand the role of the Fano effect, it is essential to take into account the coherence of 
the whole system. A model of the electron transport through a closed AB interferometer 
containing two single level quantum dots, which assumes the electrons transport through 
quantum dots are in full coherence, has been investigated in Ref. [15]. Interesting phenomena 
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such as flux-dependent level attraction and interference induced suppression of conductance 
have been found. But the effect of one-site Coulomb interaction and the interdot tunnel- 
ing haven't been considered. Ghost Fano resonance has also been observed in the study of 
electron transport through parallel DQDs with interdot tunneling but no on-site Coulomb 
interaction [3]. 

In this paper we shall investigate electron transport through parallel DQDs(schematically 
plotted in Fig.l) with interdot tunneling and on-site Coulomb interaction using the finite-?/ 
slave boson mean field theory(SBMFT) approach developed by Kotliar and Ruckenstein 17]. 
This formulation reproduces the results derived from the well known Gutzwiller variation 
wave function at zero temperature, and therefore is believed to be a powerful tool to study 
strong correlation effect of electron systems. The finite-U SBMFT has already been applied 
to investigate electron transport through single quantum do t 1 1 811 . DQDs in series in the 
Kondo regime and persistent current in mesoscopic ringj^J, and found to give good 
quantitative results for the Kondo effect on the linear conductance. 



II. THE FINITE-U SLAVE BOSON MEAN FIELD THEORY OF PARALLEL- 
COUPLED DQDS 

Electron transport through parallel DQDs with interdot tunneling and on-site Coulomb 
interation can be described by the following Anderson impurity model: 

kr\a it cr kr\ai 

(1) 

where c kr]CT {<y k ) denote annihilation(creation) operators for electrons in the leads(?7 = L,R), 
and di a (d ia ) those of the single level state in the i-th dot( i = 1,2). U is the intra-dot 
coulomb interaction between electrons, t c is the interdot tunnel coupling, and v vi is the 
tunnel matrix element between lead rj and dot i. We consider the symmetric coupling case 
with rf = rf = Yi , where Y 1 - = 2n J2 k \ v vi\ 2 ^( u ~ e fcr?o-) is the hybridization strength between 



the i-th dot and the lead r\. 

In the finite-?/ slave boson approach 







181 ] . a set of auxiliary bosons e^p^, di are intro- 



duced for each dot, which act as projection operators onto the empty, singly occupied(with 
spin up and spin down), and doubly occupied electron states on the quantum dot, respec- 
tively. The fermion operators d ia are replaced by d i<T — > fi a Zia, with z ia = e\p ia + p\ a di. 



In order to eliminate unphysical states, the following constraint conditions are imposed 
: E ct pLp«t + e\ei + d\di = 1, and fj a f ia = p\ a p i(T + djd^a =|, |). Therefore, the Hamiltonian 
(1) can be rewritten as the following effective Hamiltonian in terms of the auxiliary boson 
ei,Pi a ,di and the pesudo-fermion operators fi a : 

kr\a £=1,2 i ct 

+ X (vrri4*fl c k V * + h.c.) + x ^Q^pLput + 4 e i + d \di - 1) 

kr/ai i & 

+ ^XlHflfia-plPia-dUi), (2) 
ia 

where the constraints are incorporated by the Lagrange multipliers A 4 - and A^. The first 
constraint can be interpreted as a completeness relation of the Hilbert space in each dot, 
and the second one equates the two ways of counting the fermion occupancy of a given 
spin[17]. In the framework of the finite-U SBMFT, the slave boson operators ei,p ia ,di and 
the parameter z a are replaced by real c numbers. In this paper, we only consider the spin 
degenerate case without external magnetic field, so that all parameters are independent of 
the electron spin. We can neglect the spin index a in the parameters hereafter. Thus in the 
mean field approximation, the effective Hamiltonian is given as 

H ejf= = E^4^+ E -^fUa + t^iflf^ + H-C) 
kr)a i=l,2 & 

+ X (VmfLcbp + H - c -) + E g > (3) 

krjcri 

where t c = t c z\Z<i and v vi = v vi Zi represents the renormalized tunnel coupling between 
quantum dots and the renormalized tunnel amplitude between z-th quantum dot and the 
lead rj, respectively. Z\ and Z2 can be regarded as the wave function renormalization factors in 
the quantum dots, q = 6j + is the renormalized dot energy level and E g = £ i [A i (1) (2p? + 
e| + df — 1) — 2X^\pf + df) + Udf\ is an energy constant. 

Within this mean field effective Hamiltonian (3) the current formula through the DQDs 

m 

I = \T.\ duj[n L {uj) - n R {u)]T{u) , (4) 



is given as[ 



where the transmission probability T(u) = Tr[G a (oj)T R G r (uj)Y L ], and T L = T R 

with Ti = zfYi. The retarded/advanced Green's functions(GF) G r / a (u) 




have 2x2 matrix structures, which account for the double dot structure of the system. The 
matrix elements of the retarded GF are defined in time space as GZAt — t') = —i9(t — t') < 
{fia(t), fja-{t')} >■ By applying the equation of motion method [2 lj. one can obtain the 
retarded GF explicitly as 

-i 




zyrir 2 u - e 2 + iT 2 

The advanced GF is given by G a {u) = [G r (u)]^. Substituting the retarded/advanced GF to 
the formula of transmission probability, one obtains 



rp( , [r 1 (a;-g 2 ) + r 2 (^-g 1 ) + 2t cV /r 1 r 2 ] 2 

T \ u ) = : : : : ■ ( 6 ) 



[(a; - h){u - e 2 ) - + [T^u - e 2 ) + T 2 (lo - e x ) + 2^^^^ 

The conductance G at the absolute zero temperature in the limit of zero bias voltage is given 
byG = f \ v=o = %T(u> = 0). 

It is noticed that the formula for the transmission probability and conductance is equiv- 
alent to that of the transport through non-interacting double QD system, except that, 
in this case, the dot levels £j, the coupling strength Tj and t c are renormalized. There- 
fore, the electron transport through DQDs is characterized by the parameters e^, fj and 
t c . However, it should be noted that e^, Tj and t c show strong dependence on the gate 
voltage applied to the quantum dots, hence the result of linear conductance is quite dif- 
ferent from the non-interacting model. In the spin degenerate case, we have ten unknown 
parameters ei,pt, di, A- 1 "*, \\ 2 \i=l,2) in total to determine. From the constraints and the 
equation of motion of the slave boson operators in the effective Hamiltonian, we obtain one 
set of self-consistent equations, which is a straightforward generalized form of the single dot 
case as discussed in Ref.[18]. In this set of equations, the distribution GF of the quantum 
dots Gfj{t — t') = i < f£r{t') fur(t) > is involved, and its Fourier transform is given by 
G k (uj) = iG r {oj)\r L nL{oj) + r R riR(u})]G a (uj). We have solved the self-consistent equations 
numerically. 

In the following, we discuss the result of our calculation. First, we consider two identical 
QDs case: E\ = e 2 = and Ti = T 2 . Following Eq.(5), the transmission probability in this 
case has a Breit-Wigner resonance form, 



The retarded GF on each dot is also explicitly given as 

GI{U) = 2 1 - (e d - t c ) + 0+ + w - ( e ~ d + t c ) + 2*f ] ' (§) 

The spectral density in the i-th QD follows from the relation Piiyj) = —ImG^iu + i0 + )/n. 
It shows that the spectral density is the sum of one Lorentizan with the peak position 
located at e bond = e d + t c and one Dirac 5 peak at e antibond = e d - t c , where e bond and 
e an tibond corresponds to energy of the bonding and the antibonding state of quantum dots, 
respectively. The bonding state of DQDs has level broadening 2T due to its coupling with 
the leads. The 5 peak structure indicates that the antibonding state is totally decoupled 
from the leads. Therefore the electrons transport only through the channel of the bonding 
state, which gives a Breit-Wigner resonance form in the transmission. 

In Fig. 2 we study the effect of interdot tunneling on the transmission probability and the 
local density of state of QD in the singly occupied regime. Here we take the hybridization 
strength as the energy unit T = 1, and e d = —2. Fig. 2 (a) shows that with increasing the 
interdot coupling t c , the line shape of Breit-Wigner resonance of transmission is preserved, 
while the center of the resonance shifts to higher energy. Thus, the value of transmission 
probability at zero frequency T{uj = 0) decreases, which, in turn, results in suppression of 
the linear conductance at zero bias voltage. For local density of state shown in Fig.2(b), we 
see that the antibonding state energy is always nearby the Fermi energy of the lead, whereas 
the center of spectral density contributed from the bonding state shifts to higher energy 
along with increasing t c . 

The linear conductance G as a function of the energy level e d of the QD at zero tempera- 
ture is plotted in Fig. 3 (a) for several values of interdot tunneling t c . When there is no direct 
tunneling between two dots (t c = 0), the conductance reaches the unitary limit(G = 2e 2 /h) 
in the Kondo regime as expected. With increasing tunnel coupling t c , the conductance be- 
comes suppressed and forms a plateau structure in the regime of singly occupied QD state. 
When QDs cross over to the doubly occupied state regime, the conductance increases to the 
maximum value G = 2e 2 /h. The line shape of the linear conductance can be explained from 
the gate voltage dependence of the spectral density and the zero frequency transmission 
of QD. In Fig.3(b) we plot the bonding state energy Ibond-, the antibonding state energy 
e an tibond and the level broadening 2T as functions of e d with t c — 1. One can see that in 
the singly occupied regime with decreasing e d , the antibonding state energy e ant ibond is fixed 
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around the Fermi energy of the leads(eF = 0). It indicates that e^ ~ t c , and the conductance 
G/{2e 2 /h) = 4f 2 /[(e d + t c ) 2 + AT 2 } « l/(t 2 c /Y 2 + 1). For this identical quantum dot case, 
the value of t 2 /Y 2 is given by its bare value t 2 /Y 2 = t 2 /Y 2 . Consequently, the conductance 
shows a plateau structure and the ratio of t c /Y determines the height of the conductance 
plateau. This is in agreement with the value of conductance at the plateau structure for dif- 
ferent t c as shown in Fig. 3 (a). With a deceasing further, the QD state crosses over from the 
singly occupied to the doubly occupied regime, and €b on d goes through from positive value 
to negative value. At the point, e bond = id + t c = 0, we obtain the maximum conductance 
G = 2e 2 /h. Further down, the level broadening 2T approaches zero and the DQD will be 
totally decoupled from the leads, thus the conductance becomes zero. 

Next, we consider DQD system with different dot levels, e% ^ e 2 , and define e = (ei + e 2 )/2 
and Ae = t\ — e 2 . For the sake of simplicity, we still assume Y\ = r 2 = T. It is noted that 
in this case the renormalized hybridization strength Yi ^ Y 2 - In Fig. 4 we plot the linear 
conductance as a function of average energy of the dot levels. The parameters Ae = 0.5, 1.0 
and t c = 0.0, 0.5, 1.0, 1, 5 are used. In the case of DQDs without direct tunnel coupling t c = 0, 
Fig. 4(a) shows that the conductance curve has a narrow dip around the point e = —{7/2. 
From the formula for transmission Eq.(6), we note that the conductance vanishes when the 



condition f ie 2 + f 2 ei = 2t c y Y if 2 is satisfied. The strictly zero transmission is a consequence 
of destructive quantum interference for electron transport through the parallel DQDs, and 
it is absent for systems with DQDs connected in series. It is interesting to notice that only 
when the DQDs with different energy levels Ae ^ 0, this characteristic of interference is 
revealed. It originates from the fact that, in this case, both the bonding and antibonding 
state channels are involved in the transmission. As the energy difference Ae increases, the 
dip becomes more broadened. For non-zero interdot tunnel couplings as shown in Fig.4(b), 
(c) and (d), the conductance curves have asymmetric line shapes, which are typical for the 
Fano resonance. This results from the constructive and destructive interference processes for 
electrons transmitted through the channels of bonding and antibonding states. It is noted 
that line broadening of the Fano dip or peak depends on the value of dot level difference 
Ae, which is similar to the non-interacting DQDs caseQ]. The effect of on-site interaction 
U is to introduce strong renormalization of the dot levels and the hybridization strength, 
hence the center of the Fano resonance and line broadening have nonlinear dependence 
on the interdot tunneling t c and the level difference Ae. It is interesting to notice the 




7 



Fano resonances obtained in this study have some similarity with the experiment results in 
Ref.Q. Although their experiment is on electron transport through single QD, the coupling 
strength between quantum dot and the lead is strong and the Kondo effect and multilevel of 



22 1 have studied the 



QD might be involved in the electron transport. Recently, Biisser et al. 
electron transport through multilevel quantum dots using exact-diagonalization techniques. 
It is interesting to notice that they have also found conductance dip structure induced by 
interference effect as shown in Fig. 4(a). Actually, when t c = 0, the model studied in our 
paper is equivalent to considering two levels in single quantum dot. 

For a DQD system with energy level difference, the local density of state in i— th QD(i = 
1,2) is given by 

( , [y/r~(^ - g) + y/jy;] 2 

PiM = 1 1 —~ ■ I 9 ) 

[fa _ ~ ei ){uj - e 2 ) - i 2 ] 2 + [T x {oo - e 2 ) + T 2 (uj - e x ) + 2^^ 

In Fig. 5. we plot the local density of state in each dot. The line shape of density of state can 
be regarded as a superposition of a Fano line shape close to the antibonding state energy and 
a Breit-Wigner resonance around the bonding state energy(see Ref.[16] for detail discussion). 
The interference effect on the local density of state is manifested clearly as compared with 
that in Fig.2(b). 



III. SUMMARY 



In summary, we have studied the electron transport through DQDs in parallel config- 
uration with interdot tunneling in the Kondo regime. The strong Coulomb repulsion in 
the dots is taken into account via the finite-?/ slave boson technique. The results of our 
calculation indicate several distinct features from the non- interacting modeled : The con- 
ductance shows plateau structure as a function of the dot level in the singly occupied regime; 
Without interdot tunneling t c = 0, there is dip structure on the conductance plateau when 
the energy levels of two dots are different. When t c ^ 0, the conductance has Fano reso- 
nance line shape on the conductance plateau as a function of the averaged dot level; The 
energies of the bonding and antibonding states and the level broadening of the bonding 
state are strongly renormalized compared to the noninteracting model case. For instance, 
the antibonding state energy is almost fixed around the Fermi energy of the lead in the 
singly occupied region. The results are also different from that of the DQDs in series, in 
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which the maximum conductance is achieved when the interdot tunneling t c = 1.0 and no 



Fano resonance is observed 
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12 1. The Fano effect for parallel DQDs originates 



from the interference effect for electron transport through the two channels of bonding and 
antibonding states of parallel DQDs. In one recent experiment, Chen et al. [lOj have stud- 
ied the Kondo effect in parallel DQDs system. But the maximum conductance obtained in 
their experiment is only about 0.1e 2 /h by varying the gate voltage and interdot tunneling, 
so that we think the full coherent electron transport through DQDs isn't achieved and the 
interference effect is not manifested. One may expect further experiments on the parallel 
DQDs system can observe the conductance plateau structure and also the Fano resonance 
as discussed above. 
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Dot1 




Dot2 



FIG. 1: Parallel double quantum dots with interdot tunneling t c . 
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FIG. 2: (a) The transmission probability T{oS) and (b) the local density of state for the system 
with two identical quantum dots. Used parameters are U = 4.0, T L =T R = 1.0, e d = -2.0. The 
interdot tunnel coupling t c are 0.0(solid line), 0.5(dashed line), 1.0(dotted line). (We take the 
energy unit as T = 1, and to = corresponds to the Fermi energy of the leads). 
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FIG. 3: (a)The linear conductance as a function of the dot level at zero temperature. Used param- 
eters are U = 4.0, T L = T R = 1.0, and the interdot tunneling t c = 0.0(solid line), 0.2(dashed line), 
0.5(dotted line), 1.0 (dash-dotted line), (b) The bonding state energy e&ond (solid line), the antibond- 
ing state energy i a ntibmd (dashed line) and the level broadening 2T of the bonding state(dotted line) 
for t c = 1.0. 
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FIG. 4: The linear conductance at zero temperature for parallel double quantum dots with different 
energy levels. Used parameters are U = 4.0, T L = T R = 1.0. The energy level differences are 
Ae = 0.5 (solid line), 1.0 (dashed line), (a), (b), (c), and (d) correspond to the interdot tunneling 
t c = 0.0, 0.5, 1.0, 1.5, respectively. 
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